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Introduction: The extremely low temperatures in 
the Moon’s polar permanent shadow regions (PSR) has 
long been considered a unique factor necessary for 
entrapping volatile Hydrogen (H) [1,2]. However, re- 
cent discoveries indicate some H concentrations lie 
outside PSR, suggesting other geophysical factors may 
also influence H distributions [3-5]. In this study we 
consider insolation and its resulting thermal effects as a 
loss / redistribution process nfluencing the Moon’s 
near-surface < 1 m volatile H budget. 

To isolate regional (5° latitude band) insolation ef- 
fects we correlate two data sets collected from the on- 
going, 1.5 year long mapping mission of the Lunar 
Reconnaissance Orbiter (LRO) [6]. Epithermal neu- 
tron mapping data from the Lunar Exploration Neutron 
Detector (LEND) is registered and analyzed in the con- 
text of slope derivations from Lunar topography maps 
produced by the Lunar Observing Laser Altimeter 
(LOLA) [7][8J. Lunar epithermal neutrons are inferred 
to be direct geochemical evidence for near-surface H 
due to the correlated suppression of surface leakage 
fluxes of epithermal neutrons with increased H concen- 
tration. Regional suppressions of neutrons seen in 
LEND maps are considered localized evidence of H 
concentration increase in the upper 1m of the Lunar 
surface. To quantify spatially localized insolation ef- 
fects, LEND data are averaged from sparsely distri- 
buted pixels, classed as a function of the LOLA slope 
derivations described below. 

Background: Thermal effects of insolation are 

primarily a function of a cosine process a = i cos 0, 
which predicts the effective solar irradiation a, incident 
to a given surface as a function of its angular orienta- 
tion 0, to the source solar irradiation, i. 0 is locally a 
function of several variables including combined: sea- 
sonal, diurnal, topographic, latitude and regolith com- 
positional effects which induces locally dependent and 
time variable thermal conditions. Given the Moons 
low obliquity, increased latitude predictably attenuates 
a, and is correlated with decreased near-surface tem- 
peratures towards the poles [5], Topographic variance 
influences a slopes relative orientation to incident solar 
irradiation thus inducing local thermal conditions. Re- 
golith compositional effects e.g, granularity and com- 
position influence the thermal conductivity of the rego- 
lith near surface thus influencing the thermal gradient. 


In terrestrial polar regimes insolation effects on 
near surface distributions of water are readily con- 
trasted in pole facing (PF) vs. equator facing (EF) 
slopes. EF slopes overall receive greater expected a, 
due to their generally higher solar incident elevation 
angles vs. PF slopes. The contrasting a, thermal ef- 
fects, induce greater dessication rates tn near-surface 
EF slopes vs. equivalent PF slopes. Local PF vs. EF 
differences in a , are also maximized with increased 
slope. In this experiment we postulate similar Lunar 
effects and implement these conditions using LOLA to 
identify pixels that locally maximize a contrast and 
evaluate epithermal neutrons from LEND maps. 

Methods: LEND maps were prepared using prima- 
ry mapping mission derived LEND data (DLD) Sept 
15, 2009 to July 25, 2010. North and South polar 
maps +/-45“ to poles are produced using a 7200x7200 
pixel, 0.4 km resolution map. The mapping kernel is a 
2-D, 20km diameter uniform area mapping disk which 
is maps each sample by scaling the kernel with the 
integral of LEND’s four, 1 Hz rate, collimated sensors. 

Two polar (n,s) LOLA, 0.4 km resolution digital 
elevation models (DEM) were selected from the LOLA 
Planetary Data System release (Dec 15, 2010) files: 
ldetnj45(s,n)_400m.img, m. Image directional topo- 
graphic gradients. g K ,y = N(0,0.83)’j, y *m, where N' is 
two convolution kernels containing the (x,y) image 
coordinate directional l 51 derivatives of a 2x2 km (5x5) 
pixel. 0. centered Gaussian, o=0.83km. Topographic 
slope s = tan '(|g|). 

Slope orientation map, d> = cos"'(uv/|u||v|), defines 
a given s pixel’s angular orientation with respect to the 
pole facing direction. For each map pixel i, u, = direc- 
tional vector to the pole pixel in image space (x,v). v,- 
= slope directional vector [g x i , g y L ], The transforma- 
tion yields a linearized continuous map, range [0° to 
1 80°] of slope orientation, <J>. <S> = 0° = pole facing => 
90°=(east, west) facing => 180° = equator facing. Im- 
portantly, this is a scale invariant transformation map- 
ping all crater slopes to a common scale. This effec- 
tively facilitates geomorpholgically based classifica- 
tion. selection and sparse averaging of LEND pixels. 
Figure 1 . Due to polar trajectories, LEND map uncer- 
tainties also increase towards the equator, however 
higher low latitude uncertainties are generally offset by 
higher PF, EF averaging areas / 5° latitude bin. 




Figure 1: LOLA south pole centered DEM [8], (-80:-90). 
Elevation (greyscale). High slopes > 5 °(color). Pole facing 
f<J>:0-10) (blue). Equator facing («t>: 170-180) (red). Pixels 
for LEND/PF, EF) >= 5° slopes. 

Results: Figures 2, 3 illustrate the north, south +/- 
45° to pole comparisons of LEND pixels sparsely av- 
eraged from each 5° latitude bin. High slopes >= 5° 
are analyzed to maximize expected PF (blue) vs EF 
(red) differences in a. 5° is the upper slope threshold 
determined via pilot study for assuring each data class 
is aerially sufficiently populated. Within each latitude 
bin, LEND epithermal averages are determined as 
LOLA, fiPF = (<t> <= 5), EF = (<t> >= 170), s >= 5°). 

North results indicate epithermal neutron rate aver- 
ages collected over high EF slopes are consistently 
0.01 to 0.02 cps higher than averages collected over PF 
slopes from 50° latitude to NP. South pole results 
quantify similar EF vs. PF differences 0.01 to 0.02 cps 
for results -55° to SP. Statistical t-tests of the high 
latitude PF < EF epithermal avg. rate results indicate 
significant differences (PF vs EF) in means for a=0.05. 

Conclusion: The consistent epithermal rate (PF vs. 
EF) differences for both high latitude results +/- > 55° 
latitude suggest support for insolation effects as an 
important factor influencing dependent H loss 
processes for lunar near-surfaces < 1m in high slopes. 

We assume high slope slip faces tend to be compa- 
ratively, younger and rougher with thinner regolith 
deposits than low slope surfaces and consider the fol- 
lowing two postulates. 1 ) There exists a continuum of 
solar incidence angles between similar high PF vs. EF 
slopes. For fixed slope and roughness conditions this 
induces greater frequency and areal shadowing of small 
pockets and subsequent thermal H entrapment effects 
in PF slopes vs EF. 2) If high slope regolith is suffi- 


ciently thick, active zones similar to terrestrial perma- 
frost regimes may be produced in which volatile H 
mobility and depletion is a function of the local diffu- 
sive pentetration of heat into regolith via local a. Lo- 
cal a , heat penetration would be conditionally max- 
imized in high solar elevations near summer solstice. 

However, these are macro scale results (5° lat bins) 
and there are likely other important combinations of 
geophysical factors influencing H budgets at smaller 
scales, lower latitudes and polar contexts to consider. 

Lastly, we thank the LEND, LOLA and LRO teams 
for their instrument and spacecraft contributions to this 
and other ongoing Lunar research. 
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Figure 2(top): North Pole facing (blue: 0-10) vs. Equator 
facing slopes (red: 170-180), Latitude band average (green) 
Figure 3 (bottom): South Pole facing vs. Equator Facing 
slopes (as described in Fig 2) 



